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1. Introduction 
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It is generally believed that magnetic fields play a central role in solar 
Mh| eruptive phenomena such as flares and coronal mass ejections. The en- 

O ' ergy released through solar eruptive processes is considered to be stored 

in nonpotential magnetic fields. The magnetic energy is supplied to the 
52 I corona either by plasma flows moving around magnetic fields in the 

inertia-dominant photosphere or by magnetic flux emerging from below 
the photosphere. Since measurements of magnetic fields at the coronal 
/\ ' altitude are not available, magnetograms taken at the photospheric 

3 . level have been widely used for studies of magnetic nonpotentiality in 

flare-producing active regions and are also used through extrapolation 
to compute coronal magnetic fields. 

Several attempts have been made to identify the relationships be- 
tween time variation of nonpotentiality parameters and development 
of solar flares (Hagyard et al, 1984; Hagyard et al, 1990; Wang et 
al., 1996; Wang 1997). Moon et al. (1999c, Paper I) reviewed previ- 
ous studies on magnetic nonpotentiality indicators and discussed the 
problems involved in them. Specifically, they studied the evolution of 
nonpotentiality parameters in the course of an X-class flare of AR 6919 
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2 YONG-JAE MOON ET AL. 

using MSO (Mees Solar Observatory) magnetograms. They showed that 
the magnetic shear obtained from the vector magnetograms increased 
just before the flare and then decreased after it, at least near the 5 
spot region. Moon et al. (1999a) proposed a measure of magnetic field 
discontinuity, MAD, defined as Maximum Angular Difference between 
two adjacent field vectors, as a flare activity indicator. They applied 
this concept to three magnetograms of AR 6919 and found that the 
high MAD regions well match the soft X-ray bright points observed by 
Yohkoh. It was also found that the MAD values increased just before 
an X-class flare and then decreased after it. This paper constitutes 
one of the series of studies on evolution of magnetic nonpotentiality 
associated with major X-ray flares, which are performed using MSO 
vector magnetograms. 

Metcalf et al. (1995) studied MSO magnetograms of AR 7216, which 
are obtained from the observation of Na I 5896 spectral line, employ- 
ing a weak field derivative method (Jefferies and Mickey, 1991) and 
concluded that the magnetic field of AR 7216 is far from force-free 
at the photospheric level. This method could underestimate transverse 
field strengths for strong field regions due to the saturation effect and 
calibration problems (Hagyard and Kineke, 1995; Moon, Park, and Yun, 
1999b). Metcaff et al. (1991) used polarization data of Fe I 6302 line 
to compare two calibration methods: the weak field derivative method 
(Jefferies and Mickey, 1991) and the nonlinear least square method 
(Skumanich and Lites, 1987). They found that there are noticeable 
differences in magnetic field strength larger than B > 1200 G between 
the two methods. 

On the other hand, Pevtsov et al. (1997) analyzed 655 photospheric 
magnetograms of 140 active regions to examine the spatial variation of 
force-free coefficients. In their results, some of the active regions show 
a good correlation between B^ and J^, but others do not. It is very 
natural that the force-free coefficient varies depending on the active 
region in question even when the coefficient is more or less constant 
over the active region. Now we raise the question of whether a certain 
relation can be drawn between the force-free coefficient and the evolu- 
tionary stage of an active region. Among the selections by Pevtsov et 
al. (1997), AR 5747 is exemplary of showing a good correlation between 
Bz and Jz. Thus, we have taken AR 5747 as the object of our study. 

The purpose of this paper is to examine the magnetic nonpoten- 
tiality of AR 5747 associated with solar flares and investigate the 
evolution of the linear force- free coefficient. For this study, we have 
used the magnetograms spanning three days obtained from full Stokes 
polarization profiles of MSO. In Section 2, an description is given of the 
observation and analysis of the vector magnetograms. Computation of 
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EVOLUTION OF THE MAGNETIC FIELD IN AR5747 
Table I. Basic observational parameters of AR 5747. 



Data 


Date 


Time 


Scan 


Data Points 


Coord. 


AR5747 a) 


20 Oct., 1989 


17:41-18:50 


5.656" 


30x30 


S26W07 


AR5747 b) 


21 Oct., 1989 


19:20-20:16 


5.656" 


30x30 


S26W22 


AR5747 c) 


22 Oct., 1989 


18:27-19:41 


5.656" 


30x40 


S26W33 



nonpotentiality parameters and their evolutions in relation to flaring 
activities is presented in Section 3. In Section 4, we discuss the evolution 
of the active region field as a linear force-free field. Finally, a summary 
and conclusion are given in Section 5. 



2. Observation and Analysis 



For the present work, we have selected a set of MSO magnetograms 
of AR 5747 taken on Oct. 20-22, 1989. The magnetogram data were 
obtained by the Haleakala Stokes polarimeter (Mickey, 1985) which 
provides simultaneous Stokes I,Q,U,V profiles of the Fe I 6301.5, 6302.5 
A doublet. The observations were made by a rectangular raster scan 
with a pixel spacing of 5.6" (low resolution scan) and a dispersion 
of 25 mA/pixel. Most of the analyzing procedure is well described 
in Canfield et al. (1993). To derive the magnetic field vectors from 
Stokes profiles, we have used a nonlinear least square fitting method 
(Skumanich and Lites, 1987) for fields stronger than 100 G and an 
integral method (Ronan, Mickey and Orral, 1987) for weaker fields. 
In that fitting, the Faraday rotation effect, which is one of the error 
sources for strong fields, is properly taken into account. The noise level 
in the original magnetogram is about 70 G for transverse fields and 
10 G for longitudinal fields. The basic observational parameters of the 
magnetograms used in this study are presented in Table I. 

To resolve the 180° ambiguity, we have adopted a multi-step am- 
biguity solution method by Canfield et al. (1993) (for details, see the 
Appendix of their paper). In the 3rd and 4th steps of their method, 
they have chosen the orientation of the transverse field that minimizes 
the angle between neighboring field vectors and the field divergence 
IV-BI. 
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Table II. Basic information of the major X-ray fiares observed in AR 5747 (Solar 
Geophysical Data). 

ID Date Start (UT) End Max. Coord. Optical Class X-ray Class 

M1.4 
M2.4 
M1.9 
M3.1 
Ml. 5 
M1.3 
X2.9 



Fl 


20/10/89 


21:30 


22:03 


21:34 


S26W11 


IN 


F2 


21/10/89 


01:53 


02:06 


01:55 


S28W09 


IN 


F3 


21/10/89 


06:40 


06:49 


06:43 


S27W16 


IN 


F4 


21/10/89 


23:54 


24:00 


23:56 


S28W22 




F5 


22/10/89 


11:15 


12:30 


11:21 


S27W26 


SN 


F6 


22/10/89 


15:54 


16:47 


15:58 


S28W28 


SN 


F7 


22/10/89 


17:08 


21:08 


17:57 


S27W31 


2B 



3. Evoultion of Magnetic Nonpotentiality 

In the active region AR 5747, a number of flares took place including 
a 2B/X3 flare. In Table II, we summarize some basic features of major 
X-ray flares during the observing period. 

Figure 1 shows the ambiguity resolved vector magnetograms ob- 
tained on Oct. 20 to Oct. 22, 1989. The three magnetograms have 
the same field of view. As seen in the figures, strong sheared trans- 
verse fields are concentrated near the neutral line and they form a 
global clockwise winding pattern. In Paper I, an account is given of the 
magnetic nonpotentiality parameters used in this study. The vertical 
current density is presented in Figure 2. The vertical current density 
kernels persisted, with little change of configuration, over the whole 
observing span. Wang, Xu, and Zhang (1994) and Leka et al. (1993) 
have discussed the important characteristics of these vector magnetic 
fields and vertical current densities. We tabulate the time variation of 
magnetic fluxes and total vertical currents of positive and negative signs 
in Table III. The differences between the absolute values of the positive 
and negative quantities are within a few percent. As seen in the table, 
the magnetic fluxes and total vertical currents of both signs decreased 
with time. It is observed that several small 5 sunspots (Al, A2 and A3 
in Fig. la) disappeared in Figure lb, which suggests that flaring events 
between Oct. 20 and 21 should be associated with flux cancellation. It 
is to be noted that there were no remarkable flux emergence during the 
observing period. 

Figure 3 shows the angular shear multiplied by transverse field 
strength and Figure 4 shows the shear angle multiplied by total field 
strength. As seen in the figures, strong magnetic shear is concentrated 
near the inversion line, where H^ emission patches were observed (see 
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i)Maffnetogiam of 4Ro''47 (89/10/^0 I-4I) 




cjMdgiietogrHiti cf 4F 1 (OtJ 10 ^^ lO ^ ) 



^F( 



1000 G 




Figure 1. Three MSO vector magnetograms of AR 5747 superposed on white hght 
images taken on Oct. 20-22, 1991. In all the figures, the solid lines stand for positive 
polarities and the dotted lines for negative polarities. The contour levels correspond 
to 100, 200, 400, 800 and 1600 G, respectively. The length of arrows represents the 
magnitude of the transverse field components. 
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a) Jz of AR5747 (89/10/20 17:41) 
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# 




iltlon (aro-eo) 



Figure 2. Contours of vertical current density in AR 5747. The contour levels 
correspond to 3, 6, 9, 12 and 15 mA/m^, respectively. In all the panels, the solid 
lines stand for positive vertical current density and the dotted lines for negative 
current density. 
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a) 4ngular Shear of AR5747 (89/10/20 17:41) 




Position <«rc.oc) 



I)) 4iifriilrti Slierti of \K")717 (H'3/10/21 19 201 




Hon (.ro.eo) 



-) illglll^l Slie^ii of 4K.17 17 (BH 10/22 18:27) 







Figure 3. Contours of angular shear multiplied by transverse field strength drawn 
in thick solid lines are superposed on longitudinal magnetograms. The contour levels 
are 4.0x10*, 7.0x10", f.OxfO® and 1.3x10* Gdeg, respectively. In all the panels, the 
solid lines stand for positive longitudinal polarities and the dotted lines for negative 
ones. The contour levels in the magnetograms correspond to 100, 200, 400, 800 and 
1600 G, respectively. 
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d) Sheai Angle of AR574.7 (89/10/20 17 41) 




c) Sheai 4ii^le ni 4r>7 17 ( MH 10/33 18 37) 



a=^S' 







Figure 4- Contours of shear angle multiplied by field strength drawn in thick solid 
lines are superposed on longitudinal magnetograms. The contour levels are 4.0 x 10^, 
7.0 X 10*, 1.0 X 10^ and 1.3 x lO^Gdeg, respectively. The magnetograms are the 
same as in Figure 3. 
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c2> 





Figure 5. Contours of free energy density drawn in thick solid lines are superposed 
on magnetograms. The contour levels are 5.0 x 10*, 1.0 x 10^, 1.5 x 10^, 2.0 x 10^ 
and 2.5 x 10^ erg/cm"^, respectively. The magnetograms are the same as in Figure 3. 
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)M4D fit AKd7 17 (i^•J/lO/^G 17 11) 




— 100 - 



Position (aroseo) 

1>)M4.D of ■\Rri747 (flQ/lO/ai IQ-^O) 




: IMA.!] ot 4K.>7 1- (H<)/l()/2a Ifl 27) 



^7 




Figure 6. Contours of MAD multiplied by field strength drawn in thick solid lines 
are superposed on longitudinal magnetograms. The contour levels are 3.0 x 10*, 
5.0 X 10'', 7.0 X 10", 9.0 X lO" and 1.1 x 10^ Gdeg, respectively. The magnetograms 
are the same as in Figure 3. 
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Table III. Magnetic fluxes and total vertical currents in AR 5747 for three 
different times in the observing period. Here aj denotes the standard deviation 
of the vertical current distribution. 

Data Flux(+) [Mx] Flux(-) [Mx] Y.Jt[^] Y.Jz[H aj[mk/m^] 



a) 


1.56E22 


1.51E22 


6.7E4 


6.6E4 


1.4 


b) 


1.31E22 


1.28E22 


4.2E4: 


4.2E4 


1.2 


c) 


0.92E22 


0.89E22 


3.9E4 


4.0E4 


1.2 



Table IV. Field strength weighted mean magnetic shear angle 9s, transverse 
field weighted angular shear 6a, mean free energy density p/, planar sum of 
magnetic free energy density '}2,PSy ^^'^ X^C^^D ^ \^\) f^'' ^^ 5747 for three 
different times. The values in parentheses are obtained employing the potential 
field method. 



Data 


Os 


ea 


Pf 


Ep/ 


^(MAD X |B|) 


a) 


46.8(38.4) 


71.8(56.1) 


2.7(2.0)E4 


1.4(1.0)E24 


4.7E6 


b) 


41.3(36.3) 


63.0(53.4) 


2.1(1.6)E4 


9.3(6.9)E23 


3.5E6 


c) 


36.1(32.6) 


55.7(48.3) 


1.7(1.4)E4 


5.1(4.3)E23 


2.3E6 



Fig. 2 of Wang, Xu, and Zhang 1994). The time variation of two 
weighted mean shear angles is given in Table IV. The values of two 
shear angles monotonically decreased with time. The magnetic free 
energy density is shown in Figure 5. Its evolutionary trend is quite 
similar to that of shear angles. The 2-D MAD multiplied by total field 
strength (Figure 6) also has a similar evolutionary pattern to that of the 
other nonpotentiality parameters above. We summarize the variations 
of mean free energy density, planar sum of free energy density, and 
sum of MAD multiplied by field strength in Table IV, in which the 
values obtained with the potential field method for the 180° ambiguity 
resolution are also given in parentheses for comparison. As seen in the 
table, all the nonpotentiality parameters under consideration decreased 
with time, which suggests that the active region was in a relaxation 
stage during the observing period. 

From the above results, we may infer that the flares that occurred 
in our observation are just bursty parts of energy release in a long term 
relaxation of the stressed magnetic field. In a self-organizing system, 
a transition toward a lower energy state proceeds very mildly in the 
beginning and for most of time until a sudden bursty event develops as 
in an avalanche. Why, then, did a series of flares occur, rather than one? 
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Flares can surely take place in repetition if enough energy is supplied 
into the system between the flaring events to recover the free energy 
released by the preceding flaring event. However, this is not the case 
as far as the flares in our observation are concerned. No indication of 
energy input, whether flux emergence or increase of magnetic shear, was 
detected throughout our observing span. We thus speculate that the oc- 
currence of a series of flares was possible due to the complex geometry of 
our active region magnetic fleld. A simple bipolar magnetic fleld would 
proceed to a lower energy state by one bursty event of reconnection. 
However, in a complex active region containing more than a pair of 
magnetic poles, the transition to the lowest energy state may possibly 
comprise several steps of macroscopic change in fleld topology. This 
speculation, of course, has to be examined by further studies involving 
many other observations and numerical experiments as well. 



4. Magnetic Forces and Linear Force-Free Field 
Approximation 

It is generally believed that solar magnetic flelds are force-free in the 
corona, but far from it in the photosphere. However, to construct a 
force-free model of the coronal magnetic field, the fleld data observed 
at the photospheric level are employed as boundary conditions. Not 
only because a magnetohydrostatic equilibrium under gravity is more 
difflcult to construct than a force- free solution, but also because no reli- 
able information about plasma pressure is available in the photosphere, 
force-free fleld modeling with photospheric boundary conditions is be- 
ing widely attempted (e.g., McClymont and Mikic, 1994 for AR 5747) 
despite the afore-mentioned inconsistency. The reliability of such mod- 
els thus depends on how much the field behaves like a force-free fleld 
near the photosphere. In this section, we investigate the "force- freeness" 
of AR 5747. 

A force-free fleld is a magnetic fleld satisfying the Lorentz force-free 
condition 

(VxB)xB = 0, (1) 

which can be rewritten as 

V X B = qB . (2) 

The so called force-free coefficient a is thus given by 

Bx By Bz ' 
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in rationalized electromagnetic units. Taking divergence of Equation 
(2) and using V • B = 0, we have 

B • Va = , (4) 

which means that a is a function of each field line. With the vector 
magnetogram of Oct. 20, 1989, Canfield et al. (1991) examined whether 
the ratio of current density to field strength (J/B) is conserved along 
each elementary flux tube. Although the force- free coefficient a = Ji/Bi 
is necessarily constant along each field line in a force-free field, the con- 
dition is in practice difficult to check due to the noises in current density 
in weak field regions (McClymont, Jiao, and Mikic, 1997). To examine 
the force-freeness of AR 5747, we have computed the integrated Lorentz 
force components scaled with the integrated magnetic pressure force, 
i.e., F^/Fo, Fy/Fo and F^/Fo (Metcalf et al. 1995), in which 



1 

47r 



BxBzdxdy, (5) 



Fy = -^ I ByB.dxdy, (6) 



4tt 



Svr 
and 



j{Bl-Bl-Bl)dxdy, (7) 



^o = ^ j{Bl + Bl + Bl)dxdy, (8) 

where Fq is the integrated magnetic pressure force. In this calculation, 
only pixels with field strength larger than 100 G for both longitudinal 
and transverse fields are considered to reduce the effect of noise. In 
Table V, we present the normalized integrated forces obtained from 
three vector magnetograms. The absolute values of these forces are 
much smaller than those at the photospheric level of Metcalf et al. 
(1995), which implies that our active region field is more or less force- 
free even near the solar surface. It is also noted that the magnetic fields 
of AR 5747 become less force-free in a relaxation stage as times go. 

Now we turn to the question whether our active region field is ap- 
proximately linearly force- free. In Figure 1, the transverse field vectors 
show a common curling pattern for each magnetic polarity, which allows 
us to expect that values of the force- free coefficient do not diverge much. 
To investigate the linearity, we have plotted for each data set B^ vs. Jz 
and a plausible regression line obtained by eye fitting in Figure 7. The 
figures show that there exist approximate linear relationships between 
Bz and Jz for three vector magnetograms. We have already observed in 
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a) Bz vs. Jz of AR5747 (89/10/SO) 




t^ #*'-*+ *^ 



Is) Bz vs. Jz of AR57f7 (89/10/21) 




cQ Bz vs. .Tz of AKr)747 (»9/l0/;j:d) 




1000 2000 



Figure 7. Jz vs. Bz from three vector magnetograms of AR 5747. Straight hnes are 
the regression hnes obtained from eye fitting. 
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Table V. Normalized integrated forces and linear force free coeffi- 
cients for AR 5747. 



Data 


FjFo 


Fy/Fo 


FJF, 


Q[m-i] 


a* 


a) 


0.008 


0.005 


-0.126 


-1.1x10"^ 


-6.8x10"* 


b) 


0.042 


0.015 


-0.148 


-1.0x10"^ 


-6.1x10"** 


c) 


0.097 


0.085 


-0.222 


-7.2x10-*" 


-4.5x10"** 



Asterisked (*) values are obtained by minimizing the difference between the hori- 
zontal components of the constant a force-free model and the observed horizontal 
magnetic fields, considering only pixels with Bt > 300 G (Pevtsov et al. 1996). 



Figures 1 and 2 that the distribution of vertical electric current density 
well matches that of magnetic fluxes of opposite polarity. In Table 
V, we have tabulated linear force-free coefficients obtained by linear 
regression in Figure 7. We have also listed the coefficients obtained 
by minimizing the difference between the horizontal components in a 
constant a force-free field model and the horizontal field vectors in 
the vector magnetogram, considering only pixels with Bt > 300 G 
(Pevtsov et al, 1996). In both sets, the absolute value of force-free 
coefficients decreased with time, as other nonpotentiality parameters 
did. This suggests that the linear force-free coefficient could be as 
good a nonpotential evolutionary indicator as other nonpotentiality 
parameters as long as the linear force-free approximation is more or 
less valid. Furthermore, the linear force-free coefficient has a merit as 
a global parameter. 



5. Summary and Conclusion 

In this study, we have analyzed the MSO vector magnetograms of 
AR 5747 taken on October 20 to 22, 1989. A nonlinear least square 
method was adopted to derive the magnetic field vectors from the ob- 
served Stokes profiles and a multi-step ambiguity solution method was 
used to resolve the 180° ambiguity. From the ambiguity-resolved vector 
magnetograms, we have derived a set of physical quantities which are 
magnetic fiux, vertical current density, magnetic shear angle, angular 
shear, magnetic free energy density and MAD, a measure of magnetic 
field discontinuity. In order to examine the force- free character of the ac- 
tive region field, we have calculated the normalized integrated Lorentz 
forces and compared the longitudinal field B^ and the corresponding 
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vertical current density Jz- Most important results from this work can 
be summarized as follows. 

1) Magnetic nonpotentiality is concentrated near the inversion line, 
where flare brightenings are observed. 

2) All the physical parameters that we have considered (vertical 
current density, mean shear angle, mean angular shear, sum of free 
energy density and sum of MAD) decreased with time, which indicates 
that the active region was in a relaxation period. 

3) The X-ray flares that occurred during the observing period could 
be related with flux cancellation. Flaring events might be considered 
as bursty parts in the long term relaxation process. 

4) It is found that the active region was approximately linearly 
force-free throughout the observing span and the absolute value of 
the derived linear force-free coefficient decreased with time. Our result 
suggests that the linear force-free coefficient could be a good global 
parameter indicating the evolutionary status of an active region as long 
as the field is approximately force-free. 
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